Hard dihadron correlations in heavy-ion collisions at RHIC and LHC 
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High transverse momentum (Pt) processes are considered to be an important tool to probe 
and understand the medium produced in ultrarelativistic heavy-ion collisions via the interaction of 
hard, perturbatively produced partons with the medium. In this context, triggered hard dihadron 
correlations constitute a class of observables set between hard single inclusive hadrons (dominated 
by the leading jet fragments) and fully reconstructed jets — while they probe some features of 
the perturbative QCD evolution of a parton shower in the medium, they do not suffer from the 
problem of finding a suitable separation between soft hadrons coming from perturbative jets and 
soft hadrons coming from the non-perturbative medium as the identification of full jets does. On 
the other hand, the trigger requirement introduces non-trivial complications to the process, which 
makes the medium-modification of the correlation pattern difficult and non-intuitive to understand. 
In this work, we review the basic physics underlying triggered dihadron correlations and make a 
systematic comparison of several combinations of medium evolution and parton-medium interaction 
models with the available data from 200 AGeV Au-Au collisions at RHIC. We also discuss the 
expected results for 2.76 ATeV Pb-Pb collisions at the LHC. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

High transverse momentum {Pt) hadron production 
in the context of ultrarelativistic heavy-ion (A-A) col- 
lisions is considered to constitute an important set of 
observables to probe properties of the medium produced 
in such collisions. Often the focus is on the nuclear sup- 
pression of hard hadrons in A-A collisions compared with 
the scaled expectation from p-p collisions, which is due 
to loss of energy from the hard parton by interactions 
with the soft medium (see e.g. [l|-|5|]). The relevant ob- 
servable is the nuclear suppression factor Raa which, as 
well as its dependence on the event plane in noncentral 
collisions, has been quite thoroughly investigated in the 
context of various models [6Hll|. 

Another (visually more stunning) manifestation of the 
medium modification of high Pt partons is the appear- 
ance of hadronic mono jets in the medium if one parton 
out of a primary hard parton pair is quenched due to in- 
teractions with the medium. While such monojets have 
now been observed directly by the ATLAS and CMS col- 
laboration in 2.76 ATeV Pb-Pb collisions at the LHC 
[I2I [l3| , the phenomenon was also observed via triggered 
hard dihadron correlations by the STAR experiment at 
RHIC [l4| and more differentially measured in [l5j . 

From a theoretical perspective, the treatment of hard 
dihadron correlations is lagging behind hard single inclu- 
sive hadron production, as the correlation observable re- 
quires a substantially more involved numerical treatment 
due to the emergence of various biases which need to be 
taken into account consistently. First results [1, [TH, [l7j 



'Electronic address: thorsten.i.renk@jyu.fi 
^Electronic address: kari.eskola@phys.jyu.fi 



seemed to indicate that models which describe single 
hadron suppression well also succeed in the description of 
dihadron correlations, i.e. the results implied that there 
would not be a substantial new property of the medium 
probed in dihadron correlations. However, it was later 
realized that dihadron correlations provide tighter con- 
straints on model parameters than single hadron sup- 
pression [TU and that they are more sensitive to the 
pathlength dependence of the parton-medium interaction 
[19l ]. The dependence of the correlation on the angle of 
the back-to-back event with the bulk event has also been 
studied in 20], however no clear picture in comparison 
with data has emerged so far. 

It is the aim of this work to provide an updated view 
on the potential of hard dihadron correlations as an ob- 
servable to probe the interaction of hard partons with an 
evolving bulk medium. This is done in the light of both 
the ongoing LHC heavy ion program with a substantially 
extended reach in y/s as compared to the 200 AGeV avail- 
able at RHIC and the arrival of in-medium shower evo- 
lution Monte-Carlo (MC) codes j21H25| . Such codes are 
better suited for the treatment of the contribution of sub- 
leading jet hadrons to the correlation strength than the 
models based on leading parton energy loss which have 
often been used for the purpose so far. 



II. THE PHYSICS OF HARD DIHADRON 
CORRELATIONS 



Triggered hard dihadron correlations are interesting 
observables because they reflect an important part of 
the Quantum Chromodynamics (QCD) dynamics of a 
hard scattering process. In leading order (LO) perturba- 
tive QCD (pQCD), a hard process results in a back-to- 
back pair of highly virtual partons. Subsequently, a par- 
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tonic shower with decreasing virtuality scale develops for 
each of the two parent partons until a non-perturbative 
hadronization scale is reached and the parton showers 
turn into jets of hadrons. The total momentum balance 
of these two jets contains information about the primary 
hard QCD process whereas the momentum distribution 
of hadrons inside the jets contains information about the 
QCD dynamics of partonic shower and hadronization. 

Triggering on a hard hadron in a given momentum 
range thus on average picks the most energetic hadron of 
the two back-to-back jets. Correlations with subleading 
hadrons on the near side (i.e. in the trigger hemisphere) 
allow to probe into the QCD evolution of a jet, whereas 
correlations of the trigger hadron with the leading and 
subleading hadrons on the away side (i.e. opposite to 
the trigger hemisphere) allow to probe in addition the 
primary scattering process. 

The role of the trigger is twofold: First, it provides 
some amount of information about the kinematics at the 
hard process. In the case of 7-h or Z-h correlations where 
the QCD evolution of the trigger side is absent, the trig- 
ger momentum sets rather stringent constraints for the 
away side initial parton kinematics event by event. In 
the case of hadron-hadron correlations the constrains are 
probabilistic and the trigger momentum can only con- 
strain the parton kinematics on average. Second, select- 
ing only events in which a trigger is seen introduces a 
bias — the subclass of back-to-back hadron jets in which 
a hadron falls into a certain momentum range is different 
from the class of all back-to-back jet events. This bias 
can be systematically exploited to change the observed 
structure of the events. 

When the hard process is embedded into a hot and 
dense medium such as created in ultrarelativistic A-A 
collisions, rather general uncertainty relation arguments 
indicate that the primary hard scattering process cannot 
be substantially changed by the medium, i.e. the corre- 
sponding rate is expected to be perturbatively calculable 
and can hence serve as a 'standard candle'. However, the 
subsequent evolution of a parton shower probes time and 
distance scales characteristic for the medium evolution. 
The shower is expected to be modified by the medium 
and consequently to carry information about both prop- 
erties of the medium and the nature of parton-medium 
interaction. Finally, for light hadrons or sufficiently hard 
momentum scales, the hadronization process can be ex- 
pected to take place outside the medium. 

There is no information in hard dihadron correlation 
that cannot in principle be gained by studying full back- 
to-back jets. However, in practice in unbiased jets there is 
significant flow of energy and momentum through low Pt 
hadrons for which, in a medium, problems such as distin- 
guishing jet hadrons from medium hadrons or the possi- 
bility of in-medium hadronization (which due to its non- 
perturbative nature is poorly understood) occur. Hard 
dihadron correlations on the other hand allow to probe 
part of the QCD dynamics of jets by focusing on the 
leading jet fragments which can be required to be in a 



momentum range where such complications are absent. 



A. Observables and terminology 

The direct observables in dihadron correlation mea- 
surements are conditional yields, the so-called yield per 
trigger (YPT) on the near and away side. These depend 
crucially on the momentum range required for a trigger 
hadron. The YPT do not directly reflect any suppression 
of the observed rate of triggered events in A-A collisions 
compared to a p-p reference (which is closely related to 
the nuclear suppression factor Raa)- Thus, in an initial 
state suppression picture where a back-to-back parton 
event is either suppressed or never feels a medium, Raa 
can be at an arbitrarily low value while the YPT is un- 
changed in the medium. 

The near and away side YPT is often binned in mo- 
mentum windows, but sometimes also in terms of the 
fraction of the associate hadron momentum divided by 
trigger hadron momentum, zt = P® ssoc j P^ lg . Since zt 
has a probabilistic connection to the fractional momen- 
tum z — P^ ad /p'j? rt of a hadron produced from a parton 
with momentum p^ rt , the away side distribution binned 
in terms of zt is often called D(zt) similar to the frag- 
mentation function. Note that for 7-h correlation zt ~ z 
since the trigger photon to leading order does not frag- 
ment (it may however be produced in a shower) . 

When comparing the results from A-A and p-p colli- 
sions, usually the yield ratios Iaa — YPT me d/YPT vac 
are formed where Iaa — 1 indicates the absence of a 
medium modification to the correlation strength. Unlike 
in the case of Raa, there is no strong a priori expecta- 
tion Iaa < 1 in the medium. Various biases (discussed 
below) are potentially capable of generating Iaa > 1 for 
Raa < 1 under the right kinematical conditions. Thus, 
Iaa needs to be interpreted carefully 



B. The perturbative hard process 

In LO pQCD, the production of two hard partons k, I 
is described by 

d a AB^ki+x do l ^ kl 

, 2 , , = 2^X 1 f i/A (x 1 ,Q 2 )x 2 f 3 /B(x2,Q 2 ) £ 

dpf r dyidy 2 *-r? 1 dt 

' l J 

(1) 

where A and B stand for the colliding objects (protons 
or nuclei) and yi(2) is the rapidity of parton k(l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction x\ and a factorization scale Q ~ px is 
fi/A{x\,Q 2 ). The distribution functions are different for 
free protons [H, H3 and nucleons in nuclei |28U30| . The 
fractional momenta of the colliding partons i, j are given 
by £1,2 = ^= (exp[±yi] + cxp[±y 2 ])- Expressions for the 

pQCD subprocesses * — (s,t,u) as a function of the 
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parton Mandelstam variables s, i and u can be found e.g. 
in [p. 

To account for various effects, including higher order 
pQCD radiation, transverse motion of partons in the nu- 
cleon (nuclear) wave function and effectively also the fact 
that hadronization is not a collinear process, the distri- 
bution is commonly folded with an intrinsic transverse 
momentum kx with a Gaussian distribution, thus creat- 
ing a momentum imbalance between the two partons as 
Pt x + Pt 2 = k T . 



C. The fragmentation process 

If the full distribution of hadrons in a jet generated 
from a parent parton is known, the distribution of leading 
and up to the nth order subleading hadrons in a jet can 
be written as an expansion of the particle distribution in 
the jet in terms of a tower of conditional probability den- 
sities Ai{z\, . . . ,Zi, /i 2 ) with the probability to produce n 
hadrons with momentum fractions Z\ , . . . z n from a par- 
ton with momentum px ~ A 4 being Tl™ =1 Ai(zi, . . . Zi, ^t 2 ) 

Computing the full distributions of hadrons in a vac- 
uum or medium-modified jet is in general a compli- 
cated task requiring a model for the non-perturbative 
hadronization in addition to the perturbative evolution 
of a parton shower. In the absence of a medium, Monte- 
Carlo (MC) codes such as PYTHIA [H or HERWIG H 
are used successfully for this purpose and can be used to 
compute the Ai. 

In the medium, the Ai need to be computed from a 
MC framework for in-medium shower evolution such as 
YaJEM [22, 23] or determined approximately from a vac- 
uum shower code and a probability distribution for lead- 
ing parton energy loss P(AE). In section IIIII we will 
describe how these objects are obtained from our models 
for the parton-medium interaction in detail. Combining 
the Ai yields the full single-particle fractional momen- 
tum spectrum of produced hadrons, i.e. the medium- 
modified fragmentation function (MMFF) Djj^^z, /i 2 ). 
Here z is the fractional momentum of the hadron and n 
the momentum scale. In the leading parton energy loss 
approxination, this quantity can also be computed as the 
convolution D f ^{z, y 2 ) = P(AE) ® D f ^ h (z, p 2 ) where 
D^ h (z, /i 2 ) is the vacuum fragmentation function. 

As a side remark, while the expansion in terms of A4 
seems cumbersome when one has MC tools for the whole 
shower evolution available, it is nevertheless useful for 
two reasons. Conceptually it allows to understand the 
dependence of the subleading hadron distributions on the 
momentum of the leading shower hadron in a straight- 
forward way since e.g. Ai is a function of both Z\ and z-i. 
Technically, it allows to factorize the production of hard 
back-to-back parton pairs in the medium from the full 
in-medium shower evolution which makes the problem 
numerically easier to treat. 



D. Trigger bias in vacuum 

As mentioned above, triggering on hadrons in a certain 
momentum range is equivalent to selecting a subclass of 
events out of an unbiased distribution of back-to-back 
jets. In the absence of a medium, the trigger condition 
induces two distinct types of biases: a kinematic bias and 
a bias on parton type. 

Let us define the energy of the back-to-back parton pair 
in its own restframe as 2E. The momentum difference 
between pair restframe and cm. frame of the p-p (A-A) 
collision is then kx , leading to parton energies E\ and E2 
in the cm. frame. Assuming for simplicity for the parton 
rapidities yi = 1/2 — and negligible hadron and parton 
masses, the leading hadrons on near and away side are 
then found at Px = Z\E\ (or z\Ei\ subleading hadrons 
at z%E\ (or ZiEi). 

The kinematic bias can be understood as follows: Since 
a trigger condition Pt > Ptrig for any given primary 
parton energy E\ implies z\ > Ptrig/ Ei, for sufficiently 
small Ei of the order of Ptrig a hard fragmentation of the 
leading hadron Z\ ~ 1 is required. This in turn has impli- 
cation for the distribution of subleading shower hadrons 
Ai{z\, Z2, At) since the phase space is restricted by the 
condition Y]zj = 1. In other words, in this case the 
distribution of subleading hadrons correlated with the 
trigger softens as compared to the distribution in unbi- 
ased jets. The argument does not apply when E is large, 
however the contribution from large parton energies are 
suppressed since the pQCD parton spectrum Eq. (fT]) falls 
steeply ~ 1/E n with n — 7,8 at RHIC kinematics and 
n = 4, 5 at LHC kinematics. 

Another part of the kinematic bias concerns the mo- 
mentum imbalance cast into the form of an intrinsic kx- 
Since, in the extreme case where kx and parton momen- 
tum are parallel, for the same fragmentation function a 
parton with E\ = E + kx/2 is more likely to meet the 
trigger condition than a parton with E2 = E — kx/2, 
the trigger condition biases the distribution towards sit- 
uations in which kx is aligned with the trigger hadron 
momentum and away from situations in which kx points 
into the direction of the away side jet. The implication 
is that the momentum distribution on the away side is 
generically softer than on the near side. As in the previ- 
ous case, the bias gradually disappears for E ^ kx- 

The bias on parton type is connected with the different 
fragmentation pattern of quarks and gluons — quark jets 
exhibit a harder fragmentation pattern, i.e. are more 
likely to produce hard hadrons than gluon jets. A LO 
pQCD reaction has three possible out states: gg, qq and 
qg. Imposing a trigger condition means thus that, for 
symmetric kinematical configurations, one is more likely 
to pick events in which a quark is in the final state than 
events with a gluon. In particular, for the qg out state 
it also implies that typically the trigger hadron is part 
of the quark fragmentation whereas the away side jet 
corresponds to the gluon fragmentation. The importance 
of this latter effect depends on the Px dependent relative 
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FIG. 1: (Color online) The away side zt distribution of charged hadrons for RHIC kinematics (left panel) and LHC kinematics 
(right panel) in the absence of a medium, computed for the weight of subprocesses given by pQCD and under the assumption 
that only subprocesses with a particular out state is present. Here, the fragmentation is computed using HERWIG. 



strength of the QCD subprocess qg — > qg. 

In the following, we will distinguish between the out 
states qg and gq. We take the first case as the situation 
where the trigger hadron is part of the quark fragmenta- 
tion and the second case the situation where the trigger 
hadron was produced in the gluon fragmentation. 

The combination of kinematic and parton type bias 
has a profound (and partially counter-intuitive) effect on 
the momentum distribution of away side hadrons. Con- 
sider a kinematical situation in which the away side mo- 
mentum distribution is determined to a fraction / by 
quark fragmentation (the larger /, the harder the distri- 
bution). Naively, in a gedankenexperiment, one would 
expect that when one hardens the quark fragmentation 
the away side momentum distribution to become harder 
as well. In a certain momentum window however, harder 
quark fragmentation will imply a higher {z\) for the trig- 
ger and therefore on average smaller partonic E for the 
triggered events. However, for smaller E the bias due 
to intrinsic kr becomes more important, which tends to 
soften the away side distribution. Likewise, if the qg out 
state is a significant contribution competing with gg, then 
a harder quark fragmentation pattern will increase the 
fraction of quarks on the trigger side, implying through 
the qg — > qg processes a lower fraction f 1 of quarks on 
the away side, in essence also softening the away side 
momentum distribution. Thus, under certain kinemati- 
cal conditions, hardening the quark fragmentation may 
result in a counter-intuitive softening of the away side 
momentum distribution. 

In Fig. [1] we illustrate the combination of kinematic 
and parton type bias by showing the result for the away 
side zt distribution for RHIC and LHC for an 8-15 GeV 
trigger ran ge i n the absence of a medium, computed with 
HERWIG j34j, under the assumption that only those 



QCD subchannels which can produce the indicated par- 
tonic final state are active. For the flavour-equal out 
states qq and gg, there is very little difference in slope, 
indicating that D(zt) does not reflect the fragmenta- 
tion function D(z,fi 2 ) directly. However, for the flavour- 
unequal out states qg and gq substantial differences in 
slope are observed. At RHIC kinematics, the full pQCD 
result is chiefly an average between qq and qg which by 
accident happens to agree will with the gg result. At 
LHC, the biases are in general somewhat weaker, as the 
harder parton px spectra permit to probe larger E more 
easily. Here, the pQCD result is dominated by gg and 
softened by some admixture of qg. 



E. Trigger bias in the medium 

When a back-to-back parton pair is embedded in a hot 
and dense medium, interactions with the medium gener- 
ically soften the fragmentation pattern, and the effect 
increases with in-medium pathlength and medium den- 
sity. This has implications for the kinematic and parton 
type bias and introduces a new geometrical bias. 

Since the softening of the fragmentation increases with 
in-medium pathlength, the trigger hadron will more 
likely come from the parton which had the shorter path 
to the medium surface. In models, most partons lead- 
ing to a trigger emerge from a region ~ 4 fin from the 
medium center in central collisions [8j. This biases the 
away side parton to a long in-medium pathlength, and 
thus towards additional suppression, leading to the ex- 

, , - T away „ -i 

pectation ± AA < 1. 

On the other hand, the softening of the fragmenta- 
tion for the trigger implies (z™ ed ) < (z\ ac ) and therefore 
(E med ) > (E vac ), i.e. the mean parton energy of trig- 
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gered events is higher in the medium than in vacuum. 
Since this holds both for near and away side, for sym- 
metric flavour configurations it leads to the expectation 
Iaa V > 1 an d counteracts the geometrical bias. At the 
same time, larger (E vac ) implies that any bias due to 
intrinsic fcx is weaker than in vacuum, again leading to 

T away , 
A A > V - 

Finally, a gluon interacts by a color factor of 9/4 more 
strongly with the medium than a quark. Thus, a medium 
acts like a filter for gluons and will enhance the fraction 
of quark jets leading to a trigger hadron relative to gluon 
jets. If the out state qg is a significant contribution, this 
has implications for the away side which will then have a 
smaller fraction of quark jets than in vacuum, leading to 
the expectation I'aa V < 1 since the gluon fragmentation 
is softer. 

On the near side, smaller [z 1 [ ied ) implies more phase 
space for the production of subleading hadrons (since 

z i — 1 constrains the phase space) and hence Iaa t > 
1. The role of the medium as a gluon filter also pro- 
duces an effect into the same direction. However, the 
subleading partons in a shower are also explicitly subject 
to medium interactions which soften their distribution, 
arguing a trend I^aa? < 1 at large zt ■ 

In order to summarize and illustrate these effects in 
greater detail, let us present a simple example with semi- 
realistic numbers. Consider a hadronic trigger momen- 
tum of 10 GeV and assume for simplicity that we can 
neglect the effect of intrinsic kx- For RHIC kinemat- 
ics in vacuum, quarks are most likely to contribute to 
hard hadron production, and as a result (z\ ac ) rj 0.7, 
i.e. the typical parton momentum contributing to the 
10 GeV trigger hadron yield is ~ 14 GeV. The presence 
of a medium softens the fragmentation, and as a result 
(z™ ed ) « 0.5, thus for the same hadronic trigger momen- 
tum the typical parton momentum becomes 20 GeV due 
to the kinematic bias. At LHC kinematics, the domi- 
nant pQCD subchannel is gg — > gg and thus the frag- 
mentation pattern is softer even in vacuum. Combined 
with the harder primary spectrum of partons, one finds 
(zi ac ) 0.5, implying a parton energy of 20 GeV (con- 
siderably more than at typical RHIC conditions). The 
presence of a medium further softens the fragmentation, 
leading to (z™ 6 ^) « 0.3 and a typical parton momentum 
of ~ 33 GeV given the hadron trigger. 

The implication is that the away side before the in- 
teraction with the medium is hardened by the kinematic 
bias, and that the hardening is stronger at LHC than at 
RHIC. If one studies 8 GeV hadrons on the away side, 
one can again estimate some simple numbers using the 
parton energies estimated above. With the typical par- 
ton energies found above, one probes different regions of 
the fragmentation function, i.e. z\ ac s» 0.57 (RHIC vac- 
uum), z " led / vac ~ o.4 (RHIC medium, LHC vacuum) and 
z med py o_24 (LHC medium). For a simple exponential 
ansatz for the fragmentation function D(z) ~ exp[— 7z], 
this implies a factor ~ 9 more yield per trigger of hadrons 
at 8 GeV for LHC medium conditions than for RHIC vac- 



uum conditions before the interaction with the medium 
on the away side, a striking illustration of the kinematic 
bias connecting near side trigger bias with the away side 
momentum spectrum. The geometrical bias to have on 
average a larger in-medium pathlength on the away side 
may still be strong enough to turn this kincmatical en- 
hancement into a net suppression, but this is not a priori 
obvious. 

It is evident from these considerations that qualita- 
tive arguments are insufficient to interpret Iaa with any 
degree of confidence, as the end result will be strongly in- 
fluenced by a combination of biases acting into different 
directions. Any meaningful interpretation of experimen- 
tal data requires thus a full simulation taking all effects 
into account with their expected magnitude. 



III. THE FRAMEWORK 

The starting point for our computation of the high 
Pt hadron yield in an A-A collision is Eq. ([1]) which we 
evaluate at midrapidity y\ = y 2 = 0. We sample this 
expression using a MC code introduced in Q by first 
generating the momentum scale of the pair and then the 
(momentum-dependent) identity of the partons. A ran- 
domly chosen kx with a Gaussian distribution of width 
2.5 GeV is then added to the pair momentum. 

We assume that the distribution of vertices follows bi- 
nary collision scaling as appropriate for a LO pQCD cal- 
culation. Thus, the probability density to find a vertex 
in the transverse plane is 

P / v T A (v + h/2)T A (v -h/2) 

P{x ,yo) = ~ — 77T , (2) 

-Laa(o) 

where the thickness function is given in terms of Woods- 
Saxon distributions of the the nuclear density Pa(y, z) 
as 1a (r) — J dzpA{r,z) and TaaQi) is the standard nu- 
clear overlap function TaaO 3 ) = J d 2 sT A (s)T A {s — b) 
for impact parameter b. We place the parton pair at 
a probabilistically sampled vertex (xq, yo) sampled from 
this distribution with a random orientation <f> with re- 
spect to the reaction plane. 

We then propagate both partons on eikonal paths 
through the medium and compute for this path the lead- 
ing and first two subleading fragments using the medium- 
modified conditional probability densities Ai(zi,/j,), 
A 2 (z\,z 2 ,n) and A 3 (zi, z 2 , z 3 , m) ~ A 2 (zi + z 2 , z 3 , y) 
given the path. Finally we check if there is a hadron 
in the event which fulfills the trigger condition. If not, 
we discard the event and start generating a new one. If 
there is a trigger hadron, we bin the remaining hadrons 
in the event on the near and away side in either p^ ssoc 
or zt- 

The key point in the computation are thus the various 
medium modified Ai, which in turn require knowledge 
of the fragmentation function. For the vacuum baseline, 
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we extract the A, from either PYTHIA \M or HER- 



WIG [34| (as will be seen below, HERWIG provides a 
slightly better description of the data, but since the in- 
medium shower code YaJEM utilizes the PYSHOW algo- 
rithm [35[ , PYTHIA is the only consistent vacuum base- 
line choice if YaJEM is used for the medium modifica- 
tion). 

In the medium, the basic quantity characterizing the 
fragmentation pattern is the MMFF Dmm(z, H 2 , C) given 
the parton path. Its computation requires knowledge of 
the geometry of the medium (e.g. in terms of a space- 
time description of medium density) as well as a model 
for parton-medium interaction. Based on the results of 
we choose a 2+1 d ideal hydrodynamical model [36j 
(hydro I) and a 3+ld hydrodynamical model |37| (hydro 
II) as extreme cases for the effect of the evolving medium 
geometry for RHIC kinematics. For LHC kinematics, we 
choose an extrapolation of the 2+ld model to y/s = 2.76 
ATeV using the EKRT saturation model as described in 
11[. 

If the angle between outgoing parton and the reaction 
plane is <f>, the path of the parton through the medium 
C(r) (i.e. its trajectory £ as a function of proper medium 
evolution time r) is determined in an eikonal approx- 
imation by its initial position ro and the angle ip as 
C(t) = (xq + Tcos(0),yo + rsin(0)) where the parton is 
assumed to move with the speed of light c = 1 and the 
x-direction is chosen to be in the reaction plane. Based 
on line integrals along the trajectory, the MMFF (and 
hence the various A{) can be computed directly from the 
in-medium shower code YaJEM [H, [HI as described be- 
low. 

On the other hand, in the leading parton energy loss 
approximation the medium-modified perturbative pro- 
duction of hadrons can be computed from the convolution 

D MM (z^ 2 X) « P(AE,C)®D^ h (z,ti, 2 ) (3) 

where P(AE, Q is the medium-induced energy loss prob- 
ability given the path £ and .D* (2, p 2 ) is the vacuum 
fragmentation function. The underlying assumption is 
that the dynamics of parton-medium interactions can 
largely be cast in terms of a shift in leading parton en- 
ergy. In this case, the energy-loss probability for a given 
path £ of a parton through the medium, P(AE,Q, is 
the ingredient to be computed within a specific model of 
parton-medium interaction. 

The details of the parton-medium interaction model 
are thus in either the energy loss probability distribu- 
tion P(AE, for leading parton energy loss models or 
the MMFF Dmm{z, Mpi f° r in-medium shower models, 
given a specific path through the medium, and we will 
discuss how these are obtained in the next section. 

In all cases, the nuclear modification factor is com- 
puted given the medium-modified yield of hard hadron 
production as 



Raa{pt,v) 



T AA (b)daPP/dp T dy' 



(4) 



In the following we require for all computations that the 
parton-medium interaction model describes Raa in cen- 
tral collisions (i.e. the suppression of the rate of triggers 
in the medium) correctly both at RHIC and LHC and 
adjust model parameters accordingly. The results for the 
conditional yields are then obtained without additional 
free parameters. 



A. Armesto-Salgado- Wiedemann (ASW) formalism 

The detailed calculation of P(AE, £) follows the Baier- 
Dokshitzer-Mueller-Peigne-Schiff (BDMPS) formalism 
for radiative energy loss |38j using quenching weights as 
introduced by Salgado and Wiedemann [39j . 

In this framework, the energy loss probability 
P(AE, Q for a path can be obtained by evaluating the 
line integrals along £(r) as 



Wc (r »=/ dCCm and (qL)(v ,4>)= d(q(C) 
Jo Jo 

with the relation 



q(() = K med ■ 2 • e 3/4 (£) (cosh p - sinh p cos a) (6) 

assumed between the local transport coefficient q(Q 
(specifying the quenching power of the medium) , the en- 
ergy density e and the local flow rapidity P with angle a 
between flow and parton trajectory [40l [4l| . K me d is the 
adjustible parameter in this framework. It is naturally 
expected to be O(l), but in fits to data at 200 AGeV the 
parameter takes (dependent on the precise hydrodynam- 
ical model) values ranging between 3 and 10 (the latter 
number occurs for viscous hydrodynamics where the ini- 
tial entropy density is lower than in the ideal case, see 
.10]). ~ 

Using the numerical results of [39| and the defini- 
tions above, the energy loss probability distribution 
given a parton trajectory can now be obtained as a 
function of the initial vertex and direction (ro , </>) as 
P(AE;w c (r,<t>),R(r,<j))) = P{AE,() for w c (C) and R = 
2 Wc (C) 2 /(^(C)). 



B. AdS/QCD hybrid model (AdS) 

The AdS/QCD hybrid model is a phenomenological 
model based on the ASW formalism which is described 
in 0. 

Its basic assumption is that the hard scales in the pro- 
cess can be treated perturbatively as in the standard 
pQCD radiative energy loss calculations, while the inter- 
action with the plasma which involves strong-coupling 
dynamics can be modeled based on AdS/CFT consid- 
erations for the N = 4 super- Yang-Mills (SYM) theory 
42]. Instead of random transverse kicks, the strongly 
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coupled medium exerts a force of magnitude T 2 (with T 
the medium temperature) onto a parton and its virtual 
gluon cloud, and in the AdS model this force rather than 
transverse broadening is responsible for placing virtual 
gluons on-shell and thus leading to induced radiation. In 
a static medium of constant density, this phenomenolog- 
ical approach suggests an I? dependence of the mean 
energy loss [43| . 

While the model thus utilizes the same set of quench- 
ing weights as the ASW formalism described above, in 
order to reproduce a parametric L 3 dependence Eq. ([5]) 
is replaced by 



(g£)(r o ,0) - K AdS J dCCT 4 (C), (7) 
w c (r o ,0) = K Ads J d(( 2 T\(), (8) 
with a different model parameter KacIS- 

C. YaJEM (Yet another Jet Energy-loss Model) 

The MC code YaJEM is based on the PYSHOW code 
[H which is part of PYTHIA [33|. lt simulates the evo- 
lution from a highly virtual initial parton to a shower of 
partons at lower virtuality in the presence of a medium. 
A detailed description of the model can be found in 

mm. 

The parton shower developing from a highly virtual 
initial hard parton in this model is described as a series 
of 1 — > 2 splittings a — > be where the virtuality scale 
decreases in each splitting, i.e. Q a > Qb, Qc and the en- 
ergy is shared among the partons b, c as Eb = zE a and 
E c = (1 — z)E a . The splitting probabilities for a par- 
ton a in terms of Q a ,E a are calculable in pQCD and 
the resulting shower is computed event by event in a MC 
framework. In the presence of a medium, the main as- 
sumption of YaJEM is that the parton kinematics or the 
splitting probability is modified. In the RAD (radiative 
energy loss) scenario, the relevant modification is a vir- 
tuality gain 



d(q(C) 



(9) 



through the interaction with the medium. This modifica- 
tion leads to an increase in radiation. In order to evalu- 
ate Eq. ([IJJ during the shower evolution, the momentum 
space variables of the shower evolution equations need 
to be linked with a spacetime position in the medium. 
This is done via the uncertainty relation for the average 
formation time as 



El 
Ql 



Ea 

Ql 



(10) 



and randomized splitting by splitting by sampling rfc from 
the distribution 



P(n) = exp 



n 

(n) 



(11) 



The evolution for any given parton in the shower evo- 
lution is terminated as soon as the parton reaches a min- 
imum virtuality scale Qo- The result of the partonic evo- 
lution in terms of a shower of low virtuality partons is 
then passed on to the Lund model 44| to hadronize. The 
fractional longitudinal momentum distribution of the re- 
sulting hadron distribution corresponds to the MMFF of 
the various hadron species. 

In principle, the full functional form of q(£) could de- 
termine the MMFF, which would be computationally 
very expensive as a full MC simulation would be needed 
for every possible path in the medium. However, due to 
an approximate scaling law identified in [22], it is suffi- 
cient to compute the line integral 



(12) 



in the medium to obtain the full MMFF from a YaJEM 
simulation. The scaling law implies that the MC simula- 
tion has to be run only for a finite set of paths and makes 
a numerical solution of the geometry averaging possible. 

As in the previous case, K mec i in Eq. © serves as the 
adjustable parameter of the model once Qo is chosen. 
YaJEM requires, dependent on the underlying hydrody- 
namical model, rather natural values for K me d ranging 
from 0.6 to 2. 



D. YaJEM-D 

In the default version of YaJEM, the minimum virtual- 
ity scale is fixed at Qo = 0-7 GeV. In the version YaJEM- 
D (dynamical computation of Qo) [45| . the formation 
length of the in-medium shower is forced to be within 
the medium length. This corresponds to the choice 



Qo = sfEjl 



(13) 



which depends on both in-medium pathlength L and 
shower-initiating parton energy E. This introduces an 
explicit pathlength dependence into the MMFF and im- 
proves the agreement with the experimentally observed 
split between in-plane and out of plane emission of 
high Pt hadrons in non-central collisions [4(1 ] signifi- 
cantly. Together with the stronger pathlength depen- 
dence, YaJEM-D also predicts a stron g ris e of Raa with 
Pt in angular averaged observables |lll| which agrees 
with what has been seen by ALICE [471 ] . 
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IV. RESULTS 



B. The near side baseline 



A. The away side baseline 

In Fig. [2] we show the obtained results compared with 
the baseline d-Au data in the absence of a medium for 
three different vacuum shower MC scenarios. 

It is evident from the comparison with the data that 
all models exhibit a fragmentation pattern which is 
softer than what is seen in the data, with the default 
shower HERWIG being somewhat harder than the de- 
fault PYTHIA with Lund model hadronization. This is 
not a novel observation: In the discussion of Fig. 1 in 32 1 
we have already pointed out that the fragmentation pat- 
tern of a default PYTHIA appears somewhat softer than 
a fragmentation function such as KKP which is based on 
a global fit to data. 
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FIG. 2: (Color online) The away side zt distribution for 
RHIC d-Au collisions, computed using PYTHIA with Lund 
hadronization, PYTHIA with independent fragmentation 
hadronization and HERWIG compared with the STAR data 

El. 



In Fig. |3]we show the comparison of the vacuum shower 
evolution models with the d-Au data obtained by the 
STAR collaboration. In general, the d-Au data is de- 
scribed well by all three fragmentation functions, how- 
ever the trend is reversed with respect to what was seen 
on the away side: Here PYTHIA with Lund hadroniza- 
tion (which has the softest fragmentation pattern) leads 
to the best agreement with the slope seen in the data. 
We note that this is in line with our previous discussion 
about the kinematic bias: If the fragmentation pattern 
is in general softer, the leading hadron will on average 
appear at lower (z%), thus leaving more phase space for 
the production of subleading hadrons. 

dAu, 200 AGeV 

trigger 8-15 GeV 



Q 0.1 



0.0J 





1 1 1 




• STAR dAu data 




— PYTHIA Lund 




•— • PYTHIA indep. fragmentation - 




— • HERWIG : 







0.4 
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0.8 



FIG. 3: (Color online) The near side zt distribution for 
RHIC d-Au collisions, computed using PYTHIA with Lund 
hadronization, PYTHIA with independent fragmentation 
hadronization and HERWIG compared with the STAR data 
111. 



C. The away side in a medium 



There is some evidence that the discrepancy with the 
data largely rests in the hadronization model underlying 
the shower codes, and could in principle be removed by 
suitable tuning of parameters — if one replaces the Lund 
model by an independent fragmentation model also avail- 
able in PYTHIA, the description of the data improves. 

For the purpose of this paper, we will not attempt to 
tune the hadronization models to a best description of the 
data, but consider the output of HERWIG and PYTHIA 
with independent fragmentation as reasonably close to 
the data. Note also that the discrepancy is largely ex- 
pected to cancel in a ratio like Iaa- 



We proceed by embedding the simulation into the 
background of a hydrodynamically evolving medium. In 
order to test the sensitivity of the results to the detailed 
choice of the evolution model, we make use of the fact 
that for RHIC kinematics we with have two distinct hy- 
drodynamical models available which in a previous study 
[Io| have been found to represent extreme cases in the 
underlying geometry, leading to the weakest (hydro I) 
and strongest (hydro II) dependence of high Pt observ- 
ables on medium geometry. The differences between the 
two hydrodynamical models are largely unconnected to 
the treatment of the longitudinal coordinate as scaling 
flow (hydro I) or full simulation (hydro II). Rather, high 
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Pt observables probe the space-time distribution of ther- 
malized medium density, and hence properties such as 
the initialization time, the freeze-out conditions or the 
Equation of State being used. Here, hydro I has an early 
equilibration time of To = 0.17 fm and therefore also a 
large freeze-out temperature Tp — 160 MeV. In contrast, 
the equilibration time of hydro II is tq = 0.6 fm and 
the simulation ends at Tf = 130 MeV. Consequently, 
hydro II has a much larger freeze-out hypersurface as 
compared to hydro I, which is largely responsible for the 
observed modification of high Pt probes. For a more 
detailed comparison between the different scenarios, we 
refer the reader to [To| . 

In Fig. [4] we show D(zt) for 0-5% central Au-Au colli- 
sions at RHIC along with the appropriate baseline calcu- 
lation, computed in all four models of parton-medium in- 
teraction using the two different hydrodynamical models 
for the bulk medium. As indicated earlier, the results for 
dihadron correlations in the medium are obtained with- 
out additional free parameters. 

A few general trends can be observed: In all cases, the 
findings of [lfl| are confirmed, i.e. the hydro II medium 
evolution leads to the most pronounced medium effect. 
The leading parton energy loss models tend to underpre- 
dict the data, whereas the in-medium shower evolution 
results show a trend to be above the data. 

This is more clearly seen in Fig. [5] where we form the 
ratio with the baseline calculation and compare with the 
data in terms of Iaa(zt)- The following pattern emerges: 
In addition to the dependence on the medium model 

jhydroll < jhydrol ^ jAdS < jASW < 

jYajEM-D < jYajEM is observed . This ordering re- 
flects the pathlength dependence of the parton-medium 
interaction models, i.e. is caused by the geometrical bias, 
and agrees with the results obtained in [l(| and [45j . 
Not all combinations of medium evolution and parton- 
medium interaction model agree with the data — in par- 
ticular AdS with hydro II underpredicts the data whereas 
YaJEM-D with the hydro I overpredicts the data. In line 
with the results of 45], the pathlength dependence of 
YaJEM is too weak to lead to sizeable geometrical bias 
and hence YaJEM fails to describe the data regardless of 
the assumed medium evolution model. In contrast, the 
ASW model describes the data rather well with both hy- 
drodynamical models (ASW in combination with hydro 
I corresponds to one of the scenarios presented in our 
previous work Q). 

Let us next turn to details of the zt dependence: Here, 
the leading particle energy loss models show a quali- 
tatively wrong behaviour at low zt — while the data 
turn upward, the model results show the opposite trend. 
Given that the region zt < 0.5 probes the way sublead- 
ing hadrons in the shower are affected by the medium, 
it should not be a surprise that the leading parton en- 
ergy loss approximation is insufficient for this region (as 
indeed was already observed in [8|). The full in- medium 
shower codes on the other hand show the qualitatively 
correct upward turn at low zt, however they predict 



an onset at higher zt than seen in the data. This is 
a non-trivial constraint for the in-medium shower mod- 
elling, and indicates that the mechanism to transport 
energy towards low z in the shower is more efficient in 
nature than assumed by YaJEM. Presumably, an addi- 
tional small component of elastic energy loss is capable 
of improving the description in the low zt region. Note 
that a similar failure of YaJEM has also been observed 
in 7-h correlations 48] . 

In order to probe a different density and pathlength 
distribution, we also discuss 20-40% central collisions in 
Fig. El (for which we only have results from the 2+ld 
hydrodynamical evolution available). Here, YaJEM-D 
shows the best agreement with the data. The fact that 
YaJEM-D here slightly undershoots the data whereas for 
central collisions with the same 2+ld evolution model 
the results were above the data illustrates the strong 
non-linear pathlength dependence of the model. Clearly, 
there is currently no combination of medium evolution 
and parton-medium interaction model which gives a com- 
pletely satisfactory description of all the available data, 
thus future work on the details is needed. 

Finally, we conclude the discussion of the medium- 
modified away side by presenting Iaa(Pt) (the choice of 
Pt rather than zt is motivated by the analysis strategy 
of the ALICE collaboration) for central 2.76 ATeV Pb-Pb 
collisions in Fig. [7] This complements our discussion of 
the extrapolation of the nuclear suppression factor Raa 
from RHIC to LHC done in [Til ] . Also for this computa- 
tion, for the sake of clarity no retuning of parameters is 
done, although in the case of intrinsic kr a good case for 
a higher value can be made. Note that we only have re- 
sults from a 2+ld hydrodynamical calculation [H| (the 
extrapolation of hydro I to LHC energies) available at 
this point. It is readily apparent that Iaa° > ^aa IC ■ 
While this is surprising if one thinks in terms of geomet- 
rical bias only, the finding is readily understood in terms 
of the combination of kinematical and parton type bias. 
While for the 8-15 GeV trigger range at RHIC the qg out 
state dominates, at the LHC almost exclusively gg events 
contribute. Both the lower (zi) of the trigger hadron and 
the absence of a strong asymmetric contribution from qg 
events tend to harden the away side momentum distri- 
bution with respect to RHIC conditions. 

The remaining results can be understood in terms of 
the geometrical bias probing the pathlength dependence 
of the parton-medium interaction model. For a weak 
pathlength dependence such as found in YaJEM, the 
kinematical bias dominates and Iaa ~ 1< For stronger 
pathlength dependence there is some noticeable geomet- 
rical bias, leading to Iaa ~ 0.5 with no substantial dif- 
ferences between YaJEM-D, ASW or AdS. 



D. The near side in a medium 

It should be noted for the discussion of the near side 
correlation that unlike the away side, the dominant con- 




FIG. 4: (Color online) The away side zt distribution for RHIC kinematics for 200 AGeV dAu and 0-5% central AuAu collisions, 
compared with various combinations of parton-medium interaction and bulk evolution models. Shown are leading parton energy 
loss models (left panel) and in-medium shower evolution models (right panel) for two different hydrodynamical scenarios (see 
text). 
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FIG. 5: (Color online) The away side Iaa(zt) for RHIC kinematics for 0-5% central 200 AGeV AuAu collisions, compared with 
various combinations of parton-medium interaction and bulk evolution models. Shown are leading parton energy loss models 
(left panel) and in-medium shower evolution models (right panel) for two different hydrodynamical scenarios (see text). 



tribution to the high zt part of the near side correlation is 
given by the first subleading shower hadron. Thus, while 
the baseline calculation, being based on shower codes, is 
expected to contain all relevant physics, this is no longer 
true for the in-medium result obtained within the leading 
parton energy loss approximation. 

No significant modification of the near side associate 
hadron momentum distribution has been observed or ex- 
pected for RHIC kinematics Q, and this remains true 
for any of the models investigated here here. However, a 
near side Iaa > 1 was predicted in [32| , albeit for a very 
different trigger range, and traced back to the parton 



type bias. 

In Fig. [8] we show the results for the near side associate 
hadron momentum yield (i.e. the modification of the 
subleading shower hadron momentum distribution) for 
0-5% central Pb-Pb collisions at 2.76 ATeV. All of the 
parton-medium interaction models predict an enhance- 
ment of Iaa above unity for the 8-15 GeV trigger range, 
however the expected enhancement is much stronger in 
the in-medium shower models than in the leading parton 
energy loss models. Given that in-medium shower mod- 
els explicitly include the dynamics of subleading shower 
parton interaction with the medium whereas energy loss 
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FIG. 6: (Color online) The away side Iaa(zt) for 200 AGeV 
Au-Au collisions with 20-40% centrality, computed using a 
2+ Id hydrodynamical medium description combined with 
various models for the parton- medium interaction (see text) 
compared with the STAR data [l5j . 
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FIG. 7: (Color online) The away side Iaa{Pt) for 2.76 AGeV 
Pb-Pb collisions in the 0-5% centrality class, computed using 
a 2+ Id hydrodynamical medium description combined with 
various models for the parton-medium interaction (see text). 



models only influence this distribution via the kinematic 
and parton type bias, this finding is perhaps not a sur- 
prise. 

The observed enhancement is strongest at high Pp. 
Note that events in this bin require a substantial amount 
of energy, i.e. in addition to a subleading hadron above 
10 GeV also a triggered hadron with even larger energy 
and the energy lost from the medium. Consequently, the 
kinematic bias forces the parton energy to be very high, 
thus probing a region where the pQCD parton spectrum 
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FIG. 8: (Color online) The near side Iaa(Pt) for 2.76 AGeV 
Pb-Pb collisions in the 0-5% centrality class, computed using 
a 2+ Id hydrodynamical medium description combined with 
various models for the parton-medium interaction (see text). 

flattens and where ample phase space for secondary pro- 
duction is available. 



V. DISCUSSION 

Let us reiterate the main results of the previous section 
in order to assess the potential of hard dihadron correla- 
tions as an observable. 

• The momentum spectrum of the away side corre- 
lated yield is sensitive to both medium geometry 
and to the pathlcngth dependence of the parton- 
medium interaction model. While this is also 
true for single inclusive hadron observables such as 
Raa(<P)j the nuclear suppression factor as a func- 
tion of the angle <j> with respect to the bulk mat- 
ter event plane [l0|, the pathlength distribution in 
the background geometry is very different. Thus, a 
combined analysis of Raa(4>) an d Iaa{<P) for differ- 
ent collision centralities might help to disentangle 
the combined effects of geometry in the hydrody- 
namical model and pathlength dependence inherent 
in the parton-medium interaction model. 

• The shape of the away side zt (or Pt ) distribution 
probes energy redistribution in the shower, and is 
hence able to probe beyond leading particle energy 
loss. In the present investigation, we found con- 
straints for YaJEM which were not available from 
single hadron observables. It should however be 
noted that this set of constraints might be identi- 
cal with the ones which can be obtained from 7-h 
correlations (48j . 
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• The near side distribution likewise probes energy 
redistribution in the shower, but it does so with 
a kinematic bias in the fragmentation pattern to 
observe the trigger hadron. Thus this information 
is more differential than what can be found on the 
away side. 

• In certain kinematical regions, the parton type bias 
can be exploited to study a gluon-rich sample of 
jets. This complements the idea to study 7-h cor- 
relations in order to probe a quark-rich sample of 
jets, since qg — > 57 dominates over the annihilation 
qq 37. 

Thus, our results suggest that hard dihadron correla- 
tions are indeed a suitable tool to probe into the devel- 
opment of parton showers in the medium in a controlled 
way, starting from the leading hadron to the next sub- 
leading hadrons, while allowing for a clean separation 
between jet and medium physics at all times. 

In the present investigation, only in-medium shower 
evolution models capture correctly the qualitative be- 



haviour of the away side zt distribution seen in the data. 
Among the two in-medium shower evolution models we 
tested, only YaJEM-D remains a viable candidate since 
the pathlcngth dependence of the medium modification 
is too weak for YaJEM, and thus the away side Iaa 
overpredicts the data by a large margin. However, even 
YaJEM-D is only viable with the large spatially extended 
medium predicted in the 3+ Id hydrodynamical model. 
Thus, the hard dihadron correlations provide strong con- 
straints both for the in-medium shower models and for 
hydrodynamical modelling. 
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